Microscopic evidence for strong periodic lattice distortion in two-dimensional 

charge-density wave systems 



Jixia Dai, 1 '* Eduardo Calleja, 1 Jacob Alldredge, 1 Xiangde Zhu, 2 Lijun Li, 3 

Wenjian Lu, 3 Yuping Sun, 3 ' 2 Thomas Wolf, 4 Helmuth Berger, 5 and Kyle McElroy 1 ' ^ 

1 Department of Physics, University of Colorado at Boulder, Boulder, CO 80309, USA 
2 High Magnetic Field Laboratory, CAS, Hefei 230031, China 
3 Key Laboratory of Materials Physics, Institute of Solid State Physics, CAS, Hefei 230031, China 
4 Karlsruher Institut fur Technologic, Institut fur Festkorperphysik, P.O.B. 3640, D-76021 Karlsruhe, Germany 
5 Institute of Physics of Complex Matter, Ecole Polytechnique Federal de Lausanne, CH-1015 Lausanne, Switzerland 

In the quasi-2D electron systems of the layered transition metal dichalcogenides (TMD) there is 
still a controversy about the nature of the transitions to charge-density wave (CDW) phases, i.e. 
whether they are described by a Peierls-type mechanism or by a lattice-driven model. By performing 
scanning tunneling microscopy (STM) experiments on the canonical TMD-CDW systems, we have 
imaged the electronic modulation and the lattice distortion separately in 2H-TaS2, TaSe2, and 
NbSe2. Across the three materials, we found dominant lattice contributions instead of the electronic 
modulation expected from Peierls transitions, in contrast to the CDW states that show the hallmark 
of contrast inversion between filled and empty states. Our results imply that the periodic lattice 
distortion (PLD) plays a vital role in the formation of CDW phases in the TMDs and illustrate the 
importance of taking into account the more complicated lattice degree of freedom when studying 
correlated electron systems. 



PACS numbers: 71.45.Lr, 71.27. +a, 74.55,+v 

The transition metal dichalcogenides (TMD) have 
shown intriguing phenomena like the spin- valley physics 
[1, 2], superconductivity induced/enhanced by tun- 
ing various parameters [3-9], and charge-density waves 
(CDW) that coexist and compete with superconductiv- 
ity (e.g. 2H-NbSe 2 , TaSe 2 , lT-TaS 2 , and TiSe 2 ) [5-9]. 
Such competition between CDW and superconductivity 
is closely tied to the complicated interactions between 
the internal degrees of freedom, including charge, lattice 
and orbital. Differentiating which interactions are key 
requires knowing which of them are responsible for the 
different phases. The degree to which an electron-driven 
mechanism is the cause of CDWs in the trigonal prismatic 
structured (2H) TMDs is still under debate. Existing 
studies have proposed a variety of different mechanisms, 
including the Fermi surface nesting [10, 11], saddle band 
driven susceptibility divergence [12], f-wave gapping and 
marginal Fermi liquid [13], etc. 

The existence of incommensurate CDW phases [14, 15] 
does lend support for the role of electronic interactions as 
the cause, since 2kp is generally not expected to be a ra- 
tional fraction of the lattice reciprocal vectors. Further- 
more, the electronic origin is also supported by photoe- 
mission experiments which measured the Fermi surface of 
NbSe 2 and TaSe 2 and approximated the electronic sus- 
ceptibility with an autocorrelation. This analysis showed 
peaks at wave vectors corresponding to those of the CDW 
[10, 11]. However, inelastic X-ray scattering experiment 
in NbSe 2 [16] reveals that the lattice dynamics exhibit 
unconventional behavior and hence may be dominating 
the transition. Moreover, recent LDA calculations have 
shown the difficulties with Fermi surface nesting and have 



suggested that periodic lattice distortion (PLD), instead 
of eMod, is the essential ingredient [17, 18]. In a re- 
cent real space study, Soumyanarayanan et al. [19] have 
shown that there is a close relationship between CDW 
formation and the local strain related to the surface layer 
buckling in NbSe 2 . Rossnagel et al. have, in addition, 
argued that the spin-orbit interaction is also important 
for a 5d-electron system like TaSe 2 [20]. Because of this 
debate, it is necessary to systematically study the CDWs 
in the TMDs with the simultaneous measurement of any 
intrinsic lattice distortion while spatially resolving the 
electronic component across this family of compounds. 

The CDW instability in low dimensional systems was 
initially explained by the Peierls mechanism [21, 22] in 
which a divergence in electronic response function at a 
particular wave vector (2/cp for ID) results in a periodic 
charge redistribution in real space. That is, a low dimen- 
sional electronic system can lower its energy by opening 
an energy gap at the Fermi level, piling up electrons pe- 
riodically and reducing its lattice transitional symmetry. 
As a result of the periodic electron modulation (eMod) 
and the electron-phonon interaction, the lattice itself un- 
dergoes a similar modulation, which causes an increase 
in the lattice elastic energy [22]. As a consequence, the 
lattice distortion will be minimized and acts as a sec- 
ondary effect [18]. The eMod manifests in the change of 
the local density of states (LDOS) shows the periodic in- 
crease in density of electron with holes modulating with 
the same period but out of phase with the electrons. We 
have calculated for real-space from the Peierls model [23] . 
Another way of visualizing this alternating intensity is 
by using the sum rule in quantum mechanics, where the 
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FIG. 1. (Color online) (a) Atomically resolved topographic 
image of 2H-TaS2 with CDW modulation (3x3 superstruc- 
ture). This image was obtained at 60 K and with setup con- 
dition being —100 mV and 200 pA. Inset: zoom in to show 
the atomic structure of the 3x3 unit cell with red circles 
indicating the sulfur atoms, (b) Crystal structure of trigo- 
nal prismatic (2H) TMD. The red/black atoms are chalco- 
gen/metal atoms, and the dashed line is where it cleaves, (c) 
Fourier transform of (a). Blue (red) circles indicate the pri- 
mary peaks of CDW (atomic) modulation, and bi and 62 are 
the atomic lattice reciprocal vectors. 



summation of the numbers of filled and empty states on 
every lattice site is a constant. This means that wherever 
there are excess electron states, there must be fewer hole 
states on the same site. The LDOS calculation we have 
performed also shows this conservation of total states. 

To investigate this we use Scanning Tunneling Mi- 
croscopy (STM) in which the tunneling current is propor- 
tional to the sample's integrated density of states and a 
matrix element with exponential dependence on the tip- 
sample distance, / - e~ 2KS f* V LDOS(e)de [24]. In the 
constant current imaging mode, the tunneling current 
is maintained at a fixed value with the tip-sample dis- 
tance being adjusted and recorded. For a CDW system 
that only contains eMod and no lattice distortions, the 
topographic images of the filled and empty states will 
be complementary with each other due to the nature of 
LDOS variation and uniformity of atomic heights. This 
phenomenon has been well demonstrated by STM exper- 
iments on the quasi-lD NbSe3 [25], the 2D systems of 
Sn/Pb covered Ge(lll) surfaces [26, 27], and the quasi- 
2D CaC6 with ID modulation [28] . In each of these cases 
the modulations for filled and empty states, with the 
CDW periodicity, are of opposite phases in the topog- 
raphy (contrast inversion). 

In this study, we have used single crystal samples of 
TaS2, NbSe2 and TaSe2 prepared by iodine vapor trans- 
port method. The STM experiments were performed 
using our home-built systems, including a variable- 
temperature STM (6 to 300 K) and a low-temperature 
STM (5.3 K). In both systems, the samples were cleaved 



at room temperature in ultra-high vacuum before be- 
ing loaded into the STM. Etched tungsten tips, prepared 
in situ by field emission against gold surfaces, are used 
to obtain high quality tunneling junctions together with 
spatial resolution at atomic scale. We have examined 
three 2H-TMDs, TaS 2 , TaSe 2 , and NbSe 2 . Fig. la is 
a representative constant current topographic image of 
a 2H-TaS 2 surface showing both the atomic corrugation 
and the 3x3 CDW superlattice. Fig. lb shows crys- 
tal structure with the neutral cleave plane (dashed line) 
between the two layers of chalcogen atoms. The atomic 
lattice shown in Fig. la is the resulting triangular chalco- 
gen lattice, while the evident 3x3 superlattice is due 
to the formation of the CDW. The Fourier transform 
(Fig. lc) of Fig. la shows that the atomic and the CDW 
signals are the main features of this surface. After ex- 
cluding the atomic corrugation by Fourier filtering, we 
find that the height of the top layer of atoms vary on the 
order of picometer in the z-direction due to the forma- 
tion of CDW, for all three of the 2H-TMDs. This change 
is consistent with the distortion measured by both x- 
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FIG. 2. (Color online) Topographic images showing sym- 
metrization and antisymmetrization on 2H-TaS2 at 52 K. (a) 
and (b): Topographic images with — /+50 mV and 100 pA. 
These two images were individually corrected for drift and 
aligned with sub-atomic precision, (c) and (d): The sym- 
metrized (S = (a + b)/2) and anti-symmetrized (AS = 
(a — b)/2) images of (a) and (b). The contrast in (d) is much 
smaller compared to (c). (e): The line profiles in (a-d) indi- 
cated by colors. 
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FIG. 3. (Color online) Topographic images of NbSe2 and TaSe2 at 5 K. (a/b): Topographic images of NbSe2 with — /+100 mV 
and 100 pA. (e/f): topographic images of TaSe2 with — /+50 mV and 40/30 pA. (c/g) and (d/h): Symmetrized and antisym- 
metrized images of (a/e) and (b/f). The contrast of (c/g) is much larger than (d/h), same as TaS2- 



ray diffraction [29] and neutron scattering experiments 
[14]. This agreement strongly indicates that our surface 
sensitive measurements of the CDW show similar behav- 
ior as the bulk. 

In order to probe the differences between the filled and 
empty states of the CDW phases in these materials topo- 
graphic images were taken at opposite sample bias volt- 
ages. In Figs. 2 and 3, our topographic images obtained 
on 2H-TaS2, NbSe2 and TaSe2 are shown. Figs. 2a and 
2b are images of the filled and empty states (as indicated 
by negative and positive sample biases) in TaS2 while 
Figs. 3a-b and 3e-f are images of the empty and filled 
states in NbSe2 and TaSe2- Prior to analysis, these im- 
ages are processed with a drift-correction algorithm [30] 
so that images taken at different times can be aligned on 
a subatomic length scale. The surface defects are used to 
ensure accuracy of the drift correction algorithm. For im- 
ages of TaS2, line subtraction is applied to remove a low 
frequency vibrational noise. This process, combined with 
the high thermal stability of our experiment, allows for 
an extremely precise measurement of position (with less 
than 50 pm of drift per day) and enables us to uncover 
the true structure of the CDW state in these materials 
for the first time. 

Comparing our results with those from known CDW 
systems [25-27], which include STM results on NbSe2 
[31], we find qualitative discrepancies. In our measure- 
ments, the filled/empty states topographic images of the 
TMD samples do not show the expected out-of-phase or 
contrast inversion that the Peierls model predicts. This 
contrasts with previous observations of 2D systems such 
as Pb (or Sn) coated Ge(lll) [26, 27]. Instead, the to- 



pographic images taken at +50 and -50 mV sample bias 
voltages on TaS2 are very similar to each other when one 
compares them with the naked eye. In order to rigor- 
ously check this similarity, we have calculated the nor- 
malized cross correlation coefficients (NCC) between the 
two images, NCC = ((a^- — a) (pij — b))/ ((J a ^b) > where 
ji and a are the mean and the standard deviation of the 
two images (a^- and bij). NCC is expected to be +1 for 
two identical images, and —1 for two images with per- 
fect contrast inversion. For an eMod-dominating CDW, 
where spatial inversion between filled and empty states 
is expected, such as Pb/Ge(lll), one would expect the 
NCC between the filled and empty states images to be 
close to —1. What we find when we carry out this anal- 
ysis on TaS2, NbSe2 and TaSe2 is that the NCCs are 
+0.62, +0.79 and +0.95 respectively, showing that the 
majority of the topographic images are in-phase between 
the two biases. Considering that only the CDW corruga- 
tion will be out-of-phase while the atomic corrugation is 
actually in-phase, we can remove the atomic signals by 
Fourier filtering them and keep only the CDW signal. Af- 
ter removing the atomic signals, the measured NCC are 
still +0.37, +0.67 and +0.97 for TaS 2 , NbSe 2 and TaSe 2 
respectively, in sharp contrast to the —1 expected for a 
traditional CDW. It is important to note that this result 
differs from previous STM results on NbSe2 [31] where 
spatially out-of-phase electron and hole components were 
seen. This discrepancy is due to our superior ability to 
control and correct for drift at the subatomic level which 
allows us to uncover the true nature of these materials. 

Interestingly our measured positive valued NCCs do 
not agree with the conventional pure electronic CDW 
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model [22] and instead show that there is another non- 
electronic, lattice effect that is important to the forma- 
tion of the 2D CDW in TMDs. PLD is a natural explana- 
tion for the positive NCCs since STM measurement is a 
combination of the atomic structure (via the tip-sample 
distance) and the integrated density of states. If the ionic 
cores, around which the valence electron clouds are cen- 
tered, are displaced with the CDW transition, i.e. lattice 
distortion occurs, then the fllled-state and empty-state 
topographic images will be in-phase with each other. 
This is because they only enter into the tunneling current 
through the electron-hole symmetric matrix element. In 
other words, if there is only a periodic lattice distortion 
then the fllled-state and empty-state topographic images 
will be the same, except with a possible difference in the 
amplitudes of the height profiles, which will not affect the 
NCC. This is exactly what we see in our STM images of 
the 2H-TMDs. 

By calculating the symmetrized (Figs. 2c, 3c and 3g, 
S = (a + b)/2 where a and b are the images of filled and 
empty states) and anti-symmetrized (Figs. 2d, 3d and 3h, 
AS = (a — b)/2) images, we can directly access the PLD 
and eMod parts of the data respectively [32]. Results of 
these calculations, shown in Figs. 2 and 3 using the same 
height scaling, clearly show that the symmetrized chan- 
nel has a larger z-range than that of the anti-symmetrized 
one, agreeing with our measured NCCs above. The two 
channels (S and AS) reveal the decomposed, spatially 
resolved, PLD and eMod components of the 3 x 3 su- 
perlattice for the first time. We find similar signals for 
the two other members of the 2H family that we have 
measured, NbSe2 and TaSe2- In the framework of Peierls 
instability, the lattice distortion acts as a secondary ef- 
fect which has to be minimized because it increases the 
elastic energy of the ground state. Contrasting with the 
results of the previous STM study [31], we have demon- 
strated that the PLD is playing an important, if not the 
primary, role in the CDW formation in 2H-TMDs. This 
is calling out a strong entanglement of the PLD and the 
eMod, which is not to be expected in a purely electronic 
Peierls picture. 

Our results agree with the calculations which show that 
without the lattice involvement a CDW in these systems 
cannot exist. It has long been argued that without lat- 
tice distortions electron correlations are too strong for a 
purely electronic CDW to form [33] and that there must 
be corresponding lattice distortions. Our observations 
clearly show that these lattice distortions are present and 
dominate the observed structure in the 2H-TMDs, which 
is in contrast to previous reports to the contrary [31]. In 
addition, recent theoretical calculations in this family of 
materials have shown that the observed CDW wave vec- 
tors do not correspond to the peaks in the electronic sus- 
ceptibility [17]. This calls for a more complicated model 
that includes nontrivial phonon involvement. Such com- 
plicated origins of the charge ordering transition in such 



a simple family of materials demonstrates the necessity 
of taking the lattice contribution into account in describ- 
ing the phase diagrams of everything from the simple 
systems like TMDs, to the more complex correlated elec- 
tron materials like cuprates and manganites. 

The symmetrized topographic images do not contain 
the eMod signal because the eMod will cause asymmet- 
ric changes in LDOS [34]. Instead the symmetrized to- 
pographic images contain the PLD signal since it sym- 
metrically affects the filled and empty states. It is possi- 
ble, however, that the anti-symmetrized images do con- 
tain more information than just that of the eMod. In 
this sense, the relative intensity of PLD versus eMod we 
present is only a lower bound on the lattice contribution 
to the CDW. As we mentioned earlier, the topography is 
only sensitive to the energy window which is set by tun- 
neling bias voltage, so high energy studies might reveal 
a more complete picture of CDW in 2H-TMD. However, 
even with our low bias voltage images presented here, our 
study presents more than enough evidence to show the 
violation of the Peierls picture in the 2D CDWs in these 
materials. Furthermore, we have indeed checked these 
results by measuring over a range of bias voltages form 
25 to 200 mV and have seen no inconsistencies. 

In conclusion, we have demonstrated the importance of 
a high resolution, high stability, low drift STM to the un- 
derstanding of the structure and origins of CDWs in the 
TMDs. We have shown that atomic length scale measure- 
ments are necessary in order to resolve and separate the 
structures of both the PLD and eMod in TaS2, NbSe2 and 
TaSe 2 . This shows that the PLD dominates the CDW 
transition for the entire family of 2H-TMD and that the 
eMod, which has been previously been seen as the dom- 
inating mechanism, plays a secondary role instead. Our 
results indicate that pure electron driven CDW is un- 
likely to exist in the 2H-TMDs and that the PLD plays 
an important role, which is in excellent qualitative agree- 
ment with the recent theoretical studies [18] and provides 
quantitative benchmarks to refine our understanding. 
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LOCAL DENSITY OF STATES CALCULATION 



From the BCS-like mean field theory calculation (chapter 3 of Griiner's book, Ref. 1), 
the ground state wave function is, 

1^0) = ( n < k < k ) io> (i) 
\\k\<k F j 

where 

7i,fc = U k e-W 2 a 1<k - V k eW 2 a 2 ,k, l2,k = V k e-^ 2 a 1<k + V k e^ 2 a 2 , k , (2) 



and 



The energy dispersion is, 

E k = sign(k - k F ) (e\ + A 2 ) , e k = hv F (k - fc F ), (4) 

From ground state wave function, one can derive the density of states with a single momen- 
tum state, 

k) = 1 - 2U k V k cos(2k F x + 0), (5) 

The equations above can be found in Griiner's book [1]. 

To obtain local density of states (LDOS) which is what we measure with differential 
conductance (dl/dV), we need to convert momentum to energy, 

\E\ 



dk 

p(x,E)=p(x : k(E))— = 



1 — cos(2/c F x + (f)) 
E 



E\ 2 - IAI 2 



(6) 



The second part of this LDOS is exactly the same as the DOS predicted by BCS theory, 
while the part in the bracket is an energy dependent spatial modulation term. This term 
implies that LDOS change brought by CDW transition concentrates around E F . To visualize 
this LDOS, we have plotted a simulated line-cut (Fig. SI) and several representative LDOS's 
with different phases (Fig. S2), 

In the above simulation, the parameters we used are: fiv^ — 1, fep = 1 and hence ep = 1. 
We also used that = and A = 0.05. In order to avoid divide- by-zero issues and to include 
a physically meaningful temperature broadening, we have included a small imaginary part 
for the energy, E = E + iS. For the calculation here we used S = 0.001, which is a fair 
assumption for experiments at T>4K. 



E F 

FIG. SI. LDOS evolution with the periodicity of CDW. Three clear features are: 1. There is a 
gap around Ep; 2. States around Ep are affected the most; 3. The LDOS enhancement alternates 
between the filled and empty sides. 



25 



20 



=M5 

CO 

O 
Q 



10 











2k F x = tt/2 










2k F x = 1.2 




i 


2k F x = 0.5 




J 






2k F x = 




J 


1,1 


f 


Average 












-0.4 


-0.2 







0.2 0.4 



Energy 

FIG. S2. Representative LDOS with different phases of the CDW. The average LDOS is the same 
as the DOS predicted by BCS theory. The curves are offset for clarity. 

THEORY FOR SYMMETRIZATION OF TOPOGRAPHIC IMAGES 



I = — / [f(E F -eV + e)- f(E F + e)} LDOS(E F - eV + e) Pt (E F + e)\M\ 2 de, (7) 

tl J — OO 



The complete form of the tunneling current is as follows [2] , 
Aire 

Tl J — OO 

where f(E) is the Fermi-Dirac distribution function, LDOS and p t are density of states of 
sample and the tip, and M is the tunneling matrix element connecting the tip and sample 
wave functions. By assuming the energy window of interest is small compared to the other 
energy scales in the tunneling junction such as the work function 0, the matrix element M 
(following Bardeen, Ref. 3) is M ~ e~ KS . Here s is the tip-sample separation. Furthermore, 



we can assume T = and p t ~ const, which gives us 

I oc e~ 2KS LDOS(E F - e)de (8) 
Jo 

Thus, we see that tunneling current in a STM junction is proportional to e~ 2KS and the 
integrated density of states. 

During the constant current imaging, the tip sweeps across a sample surface with the 
feedback active. The height of the tip z{x) closely follows the height of the surface h{x) 
(defined as the position of the ionic cores and s(x) = z(x) — h(x)). By inverting the above 
equation, we get 

reV 

z{x) = h{x) + Zq \nl /l s + zq In / LDOS{x ) e)de (9) 

Jo 

From the above equation, we see that tip height should follow h(x) closely. More assump- 
tion with the integrated density of states can further reduce the above equation, 

reV 

CDOS{x ) V) = / LDOS{x ) e)de = C (V) + d(V) cos ( 9cDW x) (10) 
Jo 

where Cq is a spatial constant term and C\ is a CDW related term. Corrugation terms 
(Ci) are usually small compared to Co, so we can linearize the topographic density of states 
terms, 

z{x,V) = h(x) + z \nl /l s + z lnC Q (V) + ^cdw(^) cos(g C DW^) (H) 

where ^cdw is the amplitude of corrugation caused by CDW. This is true because any atomic 
scale variation, as verified experimentally, does not change with a change in the bias voltage 
so it can be absorbed in h{x). Except for the atomic lattice component, we expect the 
topography to only include information relating to the CDW periodicity. This is supported 
by the data we have obtained. 

In a pure electron modulation CDW, where the atomic locations are kept unchanged, 
the term ^cdw(^) will be anti-symmetric with respect to V ^cdw( - ^) -^cdw(^) = 
—za(V)- This is expected from our LDOS calculation in part I. However, if there is only a 
lattice distortion with no change in the electron-hole asymmetry, this term will be symmetric, 
zcdw(-V) = ^cdw(^) = zs(Y)- Therefore, it is reasonable to include both components 
in a general case, ^cdw(^) = za(V) + 3s (V), with zcdw(-V) = — za(V) + zs(V) where 
za and z$ are the contributions from electron modulation and periodic lattice distortion 
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respectively. It can be found that 



z(x, V) + z(x, -V) 
2 



and za(V) 



z(x, V) — z(x, —V) 
2 



(12) 



Thus, we have theoretically demonstrated how it is possible to image the electron modulation 
and lattice distortion separately by calculating these two terms. Of course if the CDW 
components were large compared to the atomic position scale or if the exponential decay of 
the states were to vary spatially these approximations would be less valid. We have verified 
that both of these approximations are valid in the systems studied here. 
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